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Abstract 

We investigate the parameter space of a supersymmetric model with an extended 
U(l) gauge symmetry in which the RH sneutrino is a thermal dark matter candidate. 
In this scenario, annihilation of RH sneutrinos proceeds mainly through Higgs or Z' 
exchange. We find that sneutrinos in the mass range from 50GeV to more than 1 
TeV can be consistent with both the WMAP limit and the direct detection upper 
limits. Powerful constraints from new gauge boson searches at the LHC as well as 
from AMg are incorporated. Depending on the choice of the U(l) charge, these 
scenarios will be further probe by direct dark matter searches as well as by Higgs 
searches at the LHC. 



1 Introduction 

The minimal supersymmetric standard model (MSSM) contains two neutral weakly in- 
teracting particles that could be dark matter (DM) candidates, the neutralino and the 
sneutrino. While the neutralino has been extensively studied and remains one of the 
favorite DM candidates [l[|5] , the left-handed sneutrino faces severe problems. The sneu- 
trino coupling to the Z boson induces a cross section for elastic scattering off nuclei that 
can exceed the experimental limit by several orders of magnitude [6] . In particular latest 
bounds from CDMS fl\ or Xenon [8] cannot be satisfied even for a LH sneutrino mass 
above ITeV. Furthermore the sneutrino annihilation rate is usually too rapid to provide 
enough dark matter [6]. 

The observation of neutrino oscillations indicative of massive neutrinos gives a natural 
motivation for adding a right-handed (RH) neutrino to the SM fields. Extending the 
MSSM with RH neutrinos and their supersymmetric partners provides then an alternate 
DM candidate, the right-handed sneutrino (RHSN). The smallness of the neutrino masses 
is usually explained by introducing Majorana mass terms and making use of the see-saw 
mechanism. The natural scale for the RH neutrinos is generally around 10^^ GeV so that 
RH neutrinos are too heavy to play a direct role in physics below the TeV scale and so 
are their supersymmetric partners. Note however that the inverse see-saw mechanism 
proposes scenarios with RH Majorana neutrinos at the TeV scale [9 , 10 . We will not 



consider these scenarios. It is also possible to generate neutrino masses through Dirac 
mass terms at the expense of introducing some large hierarchy among the fermions. In 
this case the supersymmetric partners of the neutrinos are expected to be, as for other 
sfermions, at the SUSY breaking scale, i.e. around or below ITeV. In this framework the 
RH sneutrino can be the lightest supersymmetric particle (LSP). This is the scenario we 
will consider here. 



To make a RHSN LSP a viable dark matter candidate requires special conditions. The 
RH sneutrino being sterile under standard model gauge interactions cannot be brought 
into thermal equilibrium]^ Nevertheless several proposals for sneutrino dark matter have 
emerged including mixed sneutrinos 



15 19 



RHSN in models with Dirac mass terms 

DM from a 



that result from the decay of thermal equilibrium MSSM particles 
RH sneutrino condensate 20 , sneutrinos in inverse see-saw models 



extensions of the MSSM 25 - 30 
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or RHSN in 



Extending the gauge group provides another alternative 
as the RHSN can reach thermal equilibrium because it couples to new vector and/or 
scalar fields 24 ,31,32 . For example an additional U(l)' gauge symmetry provides new 



couplings of the sneutrino with the Z' as well as to new scalar fields. In this framework 
annihilation of pairs of RH sneutrino can be efficient enough to obtain Qh"^ 0.1. The 
annihilation is specially enhanced when the particle exchanged in the s-channel is near 
resonance. Furthermore the elastic scattering cross section of the RHSN is naturally 
suppressed by several orders of magnitude as compared to the MSSM sneutrino as on the 
one hand the couplings to the EW scale particles (the Z and the light Higgs) are strongly 
suppressed and on the other hand the Z' exchange is suppressed because its mass is above 
the TeV scale. This is the framework that we will study here. 

Models with extended gauge symmetries are well motivated within the context of 
superstring models 
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grand unified theories 34 , 35 or little Higgs models 36 , 37 
Furthermore in its supersymmetric version, the additional U(l)' symmetry provides an 
elegant solution to the /i problem. Indeed as in the NMSSM [38] the effective fx parameter 
is related to the vev of the singlet S responsible for the breaking of the U(l)' symmetry [391 



40 . Another interesting feature of supersymmetric models with extended U(l) symmetry 
is that a XSHuHd interaction allows for the strong first order phase transition that is 

Finally the light Higgs mass is generally above 

Indeed 



needed for electroweak baryogenesis 41 



42 



that of the MSSM light Higgs and thus above the direct limit from LEP 
in addition to MSSM-like contributions from squarks and quarks, the light Higgs mass 
receives contributions from the superpotential (as in the NMSSM) |43j as well as from 
U(l)' D terms. 

In the MSSM with U(l) extended gauge symmetries (UMSSM), the dark matter can- 
didate can either be the neutralino or the sneutrino. The neutralino was investigated 
in (44l-[46l. We rather concentrate on the RHSN LSP. The RHSN dark matter was con- 



sidered in singlet extensions of the MSSM 27,47 in hybrid infiationnary models [48j as 
well as in different set-ups 26,49 or in models with inverse seesaw 22 . Sneutrino dark 



matter in different U(l) extensions of the MSSM was also examined in 31 and in 50 



Their signatures in cosmic rays 28 ,50 or neutrino telescopes 51 were also explored. 

In l3lj sneutrinos annihilation trough Z' exchange and B' t-channnel exchange into 
neutrinos were considered while the interactions through the Higgs sector were included 
We extend those analyses by including all possible annihilation and co- 



28,50 



m 

annihilation channels and by performing a complete exploration of the parameter space 
allowing for different choices of the U'(l) symmetry within the context of an unified 
model. Furthermore we take into account recent LHC results on the Z', on supersymmetric 
particles and on the Higgs sector. Insisting on having a sneutrino LSP, to complement 
previous studies that had concentrated on the neutralino case, will lead to some constraint 
on the parameter space. In particular the parameters and Mi,M2 have to be larger 



^Note however that non thermal mechanisms can make a mostly sterile sneutrino a good dark matter 
candidate 11-14 
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than the sneutrino mass in order to avoid a higgsino or a gaugino LSP. We will generically 
consider only cases with colored sparticles above the TeV scale to easily avoid LHC limits 
on squarks. 

Within the UMSSM with a RHSN dark matter candidate, we found a variety of annihi- 
lation channels for the sneutrino, with a predominance of annihilations near a resonance. 
The main annihilation processes can be classified as 

• Annihilation near a light Higgs resonance 

• Annihilation near a heavy Higgs resonance 

• Annihilation near a Z2 resonance 

• Annihilation into W/Z pairs through Higgs exchange 

• Coannihilation with the NLSP, here the NLSP can be neutralino, chargino, charged 
slepton or any other sfermion. 

Note that annihilation into neutrinos can occur either through exchange or through 
B\ S exchange when the latter are near the mass of the LSP. The neutrino channels are 
however never the dominant ones. Finally annihilation of sneutrinos into pairs of new 
gauge bosons could also contribute, we will not consider this channel here as it is relevant 
only for sneutrino heavier than the Z', that is above the TeV scale. For such a heavy 
sneutrino to be the LSP means that all soft parameters are also above the TeV scale. 
These configurations for supersymmetric particles have a restricted parameter space. For 
example, very large values of n often lead to a light Higgs mass below the current limits. 
Furthermore when the sneutrino is heavier than the Z', the neutralino can become the 
LSP. We therefore choose rather to concentrate on the sub TeV scale sneutrinos. For each 
scenario we also include an analysis of the direct detection rate as the spin independent 
cross section of sneutrinos on nucleons turns out to pose severe constraints on a whole 
class of U(l)' models. 

This paper is organized as follows. In section 2 we describe the various sectors of 
the model as well as the main particle physics constraints. Section 3 presents the relic 
abundance and section 4 the direct detection rate. Section 5 contains the results for two 
sample U(l)' models, the U (1)^ and U{l)n models, as well as results from a global analysis 
of the parameter space. Our conclusions are contained in Section 6, while details on the 
Higgs mass matrices are presented in the Appendix. 

2 The model 

The symmetry group of the model is SU{3)c x SU{2)l x U{1)y x U{1)' . For definiteness, 
we assume that this model is derived from an Eq model, in this case the most general 
U(l)' charges of fermions are described by only one parameter, ^^g, 

Q' = cos Oe^Qx + sin 6EsQi>, (1) 

where —n/2 < < n/2 and Qy,, Qxp are the generators of two U(l) subgroups of Eq [32] . 
The charges of the MSSM left-handed fermions as well as that of the right-handed neutrino 
are defined in Table [U 
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We also assume that the exotic (s)fermions that belong to the fundamental represen- 
tation of Eq and that are needed for anomaly cancellation are above the TeV scale. We 
will neglect them in the following. The only exception is the RH neutrino/sneutrino. 
We choose the normalization of the coupling constants such that g[ = a/5/ Sgy where 
gy and g[ are the coupling constants of f/(l)y and respectively and gy = e/cw- 

For the detailed phenomenological analysis we will consider two values of the angle 
corresponding to = 7r/2, and = — .297r {tanOEg = —a/5/3), Q.^. The first 
choice corresponds to the U(l) subgroup in the breaking of — )■ 5*0(10) x U{1)^. The 
second choice is used as an example of a scenario in which the sneutrino dark matter can- 
didate differs significantly from the case of the U{1)^ model. Note that the RH sneutrino 
decouples when 9e^ = -4:271 and that for Oei^ = the symmetry cannot be broken by the 
singlet field, S, as is the case for other U(l)' models. 



Table 1: Charges of the left-handed fermions under U'{1). 
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In addition to the MSSM superfields, the model contains a new vector superfield, 5', 
and a new singlet scalar superfield, S. The superpotential is the same as in the MSSM 
with n = but has an additional term involving the singlet, 

W = Wmssm\,=o + )^SHuH, (2) 

The vev of S, (S) = Vs/\/2 breaks the U(l)' symmetry and induces a /i term 

{I = Xvs/ V2- (3) 

The physical fields of the model are those of the MSSM with in addition a gauge boson, 
Z', and its associated gaugino, B', a scalar field, S, and a neutral higgsino, S. This 
leads altogether to 6 neutralinos and to a Higgs sector with 3 CP-even scalars, 1 CP-odd 
pseudoscalar and a charged Higgs. Each sector of the model that will play a role in dark 
matter annihilation is described below. 



2.1 Gauge bosons 



The two neutral massive gauge bosons, and Z' can mix both through mass and kinetic 



mixing 



32,44 



In the following we will neglect the kinetic mixing as it is not directly 
relevant to our analysis. The electroweak and U(l)' symmetries are broken respectively 



by the vev's of the doublets, 
reads 



'Ju/\/2 = (Hu) ,Vd/ V2 = (Hd) and singlet. The mass matrix 



Ml 



(4) 
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where 



Ml = 9[\Qlvl + Qlvl + Q'lvl) (5) 

and we have defined Qi = Q'h^i Q2 = Q'hu^ S'l ~ e/swcw- Invariance of the superpo- 
tential imposes the condition, Q'g = —{Qi + Q2)- The mixing term 

/^l = -\gig[{Qivl~Q,vl) (6) 

Diagonahsation of the mass matrix leads to two eigenstates 

Zi = cos azZ^ + sin azZ' 

Z2 = — sinazZ^ + cos azZ' (7) 



where the mixing angle is defined as 

sin2a;z 

and masses of the physical fields 



2A| 



Ml - Ml 



Ml,z, = \ (m|o + M|, T ^1 {Ml + Ml)' + A/\%^ (9) 

. The mixing angle is constrained from precise measurements of Z^ properties to be of the 
order or smaller than 10^'^ |52|, the new gauge boson Z2 will therefore have approximately 
the same properties as the Z' . As input parameters we choose the physical masses, 
Mz^ = 91.187GeV, and the mixing angle, az- From these together with the coupling 
constants, we extract both the value of tan/3 = Vu/vd and the value of Vg- Note that the 
W mass is a parameter derived from the Z^ boson mass, M\y = MzqCosOw When 
Mz2 3> Mzi and az « 1, the value of tan/3 is the solution of the following equation 

cos^/3(-g, + g,tan2/3) ^^^^sina^^ (10) 

For each U(l)' model, the value of tan /3 is strongly constrained from the lower bound 
on the Z2 mass and on its mixing. For example for U{1)^ with sina^ > 0, the value of 
tan /3 has to be below 1, this is because A| ~ 1 — tan^ /3. Large values of tan /3 are found 
when sin az < 0. The large values of tan /3 are constrained as they this might require 
too large a mixing angle for a given Z2 mass. For ^7(1),,, A| ~ 4tan^/3 — 1, therefore 
tan P < 0.5 when sina^ < 0. One might think that small values of tan /3 are problematic 
for the Higgs mass, however as we will see below additional corrections to the light Higgs 
mass can bring it above 114 GeV. 
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2.2 Higgs sector 

The CP even mass matrix is obtained after minimization of the potential. The Higgs 
sector of the UMSSM is described in l53] . We expand the Higgs field as 



H^d = ^{vd + (^d + iVd) (11) 

^n = ^ i^u + 0n + ^^u) (12) 

S=^{v, + a + tO, (13) 

The squared mass in the (j)d,4>u,<^ basis both at tree level and including the dominant 



radiative corrections from top quarks and stops is listed in Appendix A 53 . This 
squared mass matrix is diagonalized by the transformation Md = ZhAiZj^^ and the mass 
eigenstates (ordered in mass) are {hi,h2,h^Y = Z'[{(f)d,4>u,o')'^. 

There is only one pseudoscalar Higgs which is obtained after diagonalizing the 3x3 
pseudoscalar mass matrix listed in Appendix A. This matrix is diagonalized by the unitary 
transformation Mjj = ZaM-a^a^ with the physical pseudoscalar, = ZAi^Pd + ZA2fu + 
ZA3i- The tree- level mass 

rriA = -V 14 

^ sin 20 ^ ^ 

and tan0 = t>sin2/3/2t;s. When Vg » v, the pseudoscalar mass reduces to 

tti a ~ - — (15) 

^ sm/3cos/3 ^ ^ 

The full matrix including one-loop corrections is also listed in Appendix A together with 
the corrections to the charged Higgs mass. 

Typically the Higgs spectrum will consist of a standard model like light CP-even Higgs, 
a heavy mostly doublet scalar which is almost degenerate with the pseudoscalar and the 
charged Higgs, and a predominantly singlet scalar. Note that the mass of the pure singlet 
~ g'lQsVs and is therefore close to that of Mz^ when Vg 3> Vu, Vd, see Eq. [4j The 
hierarchy in the mass of the heavy doublet and the singlet depends on the parameters 
of the model. For large values of A\ and A (and therefore fi) the mass of the heavy 
Higgs doublet increases and can exceed that of the singlet. Furthermore the Higgs mixing 
increases, however the singlet component of the light state is usually not large enough to 
significantly relax the bound on the lightest Higgs contrary to what occurs in the NMSSM 



38 ,54 56 . Although the lightest scalar Higgs is usually standard model like, it can be 



significantly heavier than in the MSSM. Indeed the upper bound 53 on the Higgs receives 
two types of additional contributions as compared to the MSSM, one proportional to Ax 
that is also found in the NMSSM, and one from the additional gauge coupling g[. The 



upper bound on the lightest Higgs mass is thus raised to (9(170) GeV 53 . Such a large 
Higgs mass is however constrained by the latest LHC results which for standard model 
couplings exclude ruh > 144 GeV at the 90%C.L. [571. The impact of the LHC results on 



the parameter space of the model will be discussed in more details in 58 . 
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2.3 Sfermions 



The important new feature in the sfermion sector is that the U(l)' symmetry induces 
some new D-terms contributions to the sfermion mass, Af, 



where Q'j are the U(l)' charges of the left fields listed in Table [TJ The sfermion mass 
matrix reads 



tribution can completely dominate the sfermion mass, especially for large values of Vs- 
Those are found in particular when ~ 0. The D-term contribution can induce neg- 
ative corrections to the mass, so that light sfermions can be found even when the soft 
masses are set to 2TeV. For ^£;g > 0, Aj, < Ag^uR^en < ^L,df, so that a universal soft mass 
term for the sfermions at the weak scale naturally leads to a RH sneutrino as the light- 
est sfermion. Furthermore the NLSP will be the right-handed slepton or the t-squarks 
if a large mixing decreases the mass of the lightest t. On the other hand for Oe^^ < 0, 
Au > Aq^uRfin > ^L,dii SO the sneutrino cannot be the LSP with universal soft sfermion 
masses at the weak scale. In general one expects non-universality in sfermion masses at 
the weak scale even if universality is imposed when the model is embedded in a GUT 
scale model. In particular the right-handed sleptons, whose RGB's are driven only by 
U(l) couplings have the smallest soft terms at the weak scale. Therefore as long as Ay 
is not much larger than for other sfermions, it is still natural that the sneutrino be the 
lightest sfermion. This is the case for models where Oe^ ^ —tt/2 since all correction terms 
are then almost equal. For the sake of reducing the number of free parameters we will in 
this study assume that the soft sfermion masses are universal at the weak scale safe for 
that of the RH sneutrino. In some cases we will find a sfermion NLSP that contributes 
to coannihilation processes. In the U{l)rf model the NLSP will be mostly the left-handed 
slepton or down-type squarks while in the U{1)^ model it will be rather right-handed 
slepton or up-type squarks. 

2.4 Neutralinos 

The mass terms for the gaugino sector reads 



A/ = -g?Q'f {v\Q, cos' /3 + g2 sin^ /3) + Q'^v',) 



(16) 



2 I 




) 



where Qf and 




(18) 
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The neutralino mass matrix in the basis {B, W^, H^, H^, S, B') thus reads 

^ Ml -Mzoc^sw Mzos/ssw Mk ^ 

M2 Mzoc,3Cw -Mzospcw 



— MzoCpSw MzoCjjCw 

MzoSisSw -MzoSpCw 



Mk 



-/i -^ivsp/vs Qig'ivc/s 

-/i -flVCp/Vs Q2g'iVSfi 

-fiVS^/Vs -fiVCfj/Vs Q'sg'lVs 

Qxg'xvcfi Q2g'\VSfi Q'sg'iVs M[ J 



The mixing matrix is defined such that the hghtest neutrahno is written as 
X? = NuB + N,,m + m^Hd + N^aHu + iVi5^ + N,^B> 



(19) 



(20) 



Note that >> f so that one expects a large mixing between the singhno and B' . The 
chargino sector is identical to that of the MSSM. 



2.5 Constraints on the model 

The main constraints on the model arise from the gauge boson and Higgs sector. Since 
colored supersymmetric particles do not play a direct role in sneutrino annihilation, we will 
generally assume that they are above the TeV scale thus evading the LHC constraints 59 , 
60|. 

The Z2 can be produced directly in pp and pp collisions and is searched primarily 
using leptonic decay modes. The limits from the LHC released earlier this year have 



superseded the Tevatron hmits |52 
in 



We will use the ATLAS exclusion limits published 
61 obtained with an integrated luminosity of £ = 1.01(1.21)/6~^ in the e"'"e~(/i'^ 



/i 



channels. These limits are extracted using different U(l)' models, in the two models we will 
consider for the phenomenological analysis the bounds are similar, with Mz^ > 1.49TeV 
for U{1)^ and > 1.54 TeV for f/(l)r,- These bounds were derived assuming the Z2 
decays only into SM particles. In our case, the Z2 can also decay into supersymmetric 
particles, into RH neutrinos and into Higgses, thus reducing the branching ratio into e, /z. 
For example in the U{1)^ model the decays into supersymmetric particles can reach up 
to 20% especially in cases where the sneutrino is light (hence neutralinos can be light 
as well) while decays into neutrinos are typically below 10%. In the U{l)rj, the decays 
into SM particles is even more suppressed, it never exceeds 65%. This is mainly due to a 
large branching fraction into right-handed neutrinos ( around 30%) as well as into other 
sparticles. The limits on the Z2 mass are therefore weakened. To take this effect into 
account we have computed the modified leptonic branching ratio for each point in our 
scans and have re-derived the corresponding ATLAS limit. 

Another important constraint on the gauge boson sector comes from precision mea- 
surements at the Z-pole and from low energy neutral currents. These provide stringent 
constraints on the mixing angle az- Depending on the model parameters the constraint 
are below a few 10" 



52 



We will typically set the mixing angle to \az\ = 0.001. 
The Higgs sector is also severely constrained, a standard model like Higgs is constrained 
by LEP to be above 114 GeV [62j and by LHC searches to be below 144GeV [57j. These 
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limits can be relaxed when the Higgs couplings to gauge bosons is reduced due to the 
mixing with the singlet. In the parameter space explored, the singlet component is small 
and the limit is not modified significantly. The upper limit can also be relaxed if the Higgs 
has a large branching fractions into invisible particles. The relic density constraint imposes 
that light sneutrinos have a mass very near Despite a phase space suppression 

factor, the contribution of the invisible mode to the light Higgs decay can in some cases 
reach 90% for light sneutrinos, thus relaxing the constraint on the lightest Higgs. When 
imposing a limit on the light Higgs mass, we have folded in the effect of the invisible decay 
modes therefore allowing in a few light Higgs heavier than 144 GeV. 

Observables in the B sector provide powerful constraints on the supersymmetric pa- 
rameter space assuming minimal flavour violation. The mass differences of B mesons, 
AMg, AMrf have been measured to be 17.77 ± .12ps~^ and 0.507 ± .004ps~^ 62 , some- 
what below the SM predictions 

AMf*^ = 20.5±3.1j9s-i 

AM|*^ = 0.59±0.19ps"^ (21) 

Additional supersymmetric contributions of the same sign as the SM ones are there- 
fore strongly constrained, even though there are large uncertainties in the SM prediction 
mainly due to the CKM matrix elements and hadronic parameters. Supersymmetric con- 
tributions include box diagrams with charged Higgs and quarks, squark/chargino, neu- 



tralino and or gluino loops 63 as well as double penguin diagrams with a neutral Higgs 
exchange. The latter give a significant contribution for large values of tan /5 and were not 
included in our analysis. The former are important at small values of tan /3 which are often 
found in the U(l)' model. In particular the charged Higgs/quark box diagram contribution 
adds to the SM contribution inducing too large values for AM^. The dominant contri- 
bution, for small values of tan/3 is proportional to a;logx(cot^/3), where x = m^/mjj+. 
This observable thus constrains severely scenarios with tan /3 < 1. The computation of 



the mass difference is adapted from the routine provided in NMSSMTools 64 . We have 
also used the same estimate for the theoretical uncertainties. 

Other observables such as 6 — )■ 57, Bs — )• /i"*"//" or B ^ tu are known to receive large 
supersymmetric contributions when tan /3 is large, the heavy Higgs doublet is light and/or 
the squarks are light. The scenario we will study have heavy squarks, Higgs doublets above 
several hundred GeV's and feature small to moderate values of tan f3, we therefore have 
not included these constraints. 



3 Relic abundance of sneutrinos 

The RH sneutrino couples very weakly to the MSSM particles, its annihilation is therefore 
primarily through the particles of the extended sector, the new vector boson, Z2 and new 
scalars. The couphng of the RHSN to Z2 is directly proportional to its U(l)' charge, 

gz2VRU*^ = g'lQ'u cos az (22) 

This coupling does not depend much on the other parameters of the model since cos az ~ 
1. Annihilation of sneutrinos becomes efficient when m^^ ~ M^2/2 since the s-channel 
processes can benefit from resonance enhancement. 
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Sneutrino annihilation also become efficient when m/j^ ~ ^kI"^ where hi can be any 
neutral CP even Higgs field. The couphngs of the sneutrino to neutral scalars reads 

9h,vnvR = -g'lQ'v b cos pZhiiQi + V sin ^Zhi2Q2 + VsZhisQ'sl ■ (23) 

Since f , the largest coupling will be to the predominantly singlet Higgs, for which 

Zh3i ~ 1. Typically the singlet Higgs is the one that has a mass close to Z2, resonant Higgs 
annihilation will therefore occur for roughly the same sneutrino mass as the resonant 
annihilation through Z2. The light Higgs is dominantly doublet, nevertheless its coupling 



to sneutrinos receives contributions from all three terms in Eq. 23 since Vs ^ v. This 
coupling is generally sufficient to have a large cross section enhancement near rrij^^ ~ 
In some cases efficient annihilation can occur away from the resonance, for example 
annihilation into W/Z pairs through light Higgs exchange or - for heavier RHSN's - into 
light Higgs pairs or ti pairs through singlet exchange. Note that the sneutrino coupling 
to the lightest Higgs depends on A (/x), an important parameter to determine the mixing 
of Higgses. For some choice of parameters, the couplings of the mostly doublet hi to 
the RHSN can be strongly suppressed, not allowing a large enough enhancement on the 
annihilation cross section. Specific examples will be presented in section [5) 

Sneutrinos can also annihilate in neutrino pairs through t-channel exchange of B' in 
addition to the usual Z2 contributions. This process contributes mostly for light B' and 
is never the dominant channel. Finally it is also possible to reduce the abundance of 
sneutrinos through coannihilation processes, this occur when the masses of the NLSP and 
the LSP are within a few GeV's. Coannihilation can occur with neutralinos, charginos 
or other sfermions. Typically, because there is only weak couplings of the RH sneutrino 
to the rest of the MSSM particles, coannihilation processes involve self-annihilation of 
the NLSP's and/or NLSP/NNLSP annihilation. The neutralino NLSP will decay via 
Xi ~^ ^R^R with a typical lifetime of 10"^^ — 10^^''s except when it is almost degenerate 
with the sneutrino LSP which can lead to an increase of lifetime. The NLSP decay will 
always occur much before BBN and will not spoil its predictions. 



To compute the relic abundance we use micrOMEGAs2.4 65 . For this we first imple- 



10 



ment the model in LanHEP [66' in both unitary and Feynman gauges and checked gauge 
invariance for a large number of processes at tree-level. We then introduced radiative 
corrections to the Higgs masses in the unitary gauge. This code then produces the model 
file for CalcHEP 



67 and micrOMEGAs 68 



All facilities of micrOMEGAs for computing 
dark matter properties are then obtained automatically, in particular the relic density 
and the LSP-nucleon scattering cross section. When describing the numerical results we 
will always impose the condition that the sneutrino is the LSP. For another study of the 
neutralino LSP as a DM candidate see 



44 



4 Direct detection 




hi 




q q q q 

Figure 2: Main scattering processes of the RH sneutrino. 



The cross section for scattering of sneutrinos on nucleons will be purely spin indepen- 
dent as the sneutrino is a scalar particle. This process receives contributions from only 
two types of diagrams: exchange of gauge or Higgs bosons. The gauge boson contribution 
depends on the vectorial coupling of the fermions to Zi^2, with Qy = Is — 2Qsm^ 9w for 
Z° and Qy = Q'^ — Q'jc for Z'. The total gauge boson contribution to the spin independent 
cross section on point-like nuclei reads 



4 

TT 

where 



= —{g'lQ'.Y - 4sin2^H.) +y')Z + {-y + 2y'){A - Z)]' (24) 



and Qy = — 4/ a/40 cos 6'^;^. The contribution in l/Mf^ in y and the one proportional to 
sin^ az in y' are always suppressed. From Table [l] one sees that all fermions have a purely 
axial-vector couplings to Z^ and that u-quarks also have purely axial vector couplings to 
Z^. Therefore the Z' contribution is solely dependent on its couphng to d-quarks, hence 
the term in Qy in y'. This contribution is expected to be twice as large for neutrons than 



for protons, see Eq. 24 In the model f/(l)^, y' = and the cross section will depend on 



the Zi exchange contribution. In this case the amplitude for protons is suppressed by a 
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factor 1 — 4 sin^ 6w as compared to that for neutrons. Furthermore the Zi contribution 
proceeds through the Z' component so that the cross section is proportional to sin^a^. 
When cos6'£;g 7^ the term in y' usually dominates by about one order of magnitude for a 
TeV scale Z2 and the mixing angle az = 10^^. In these scenarios only a weak dependence 
on az is expected, the largest cross sections are expected for Oe^^ — ^ corresponding to 
the maximal value of Qy. 

The Higgs contribution leads to a cross section 

-fA = ^ E f - ^) E J>n + zY: 9..,Jl^^ , (26) 

where ghiqq = —e-ZhnMq/ {2M-^ySwCi3) is the Higgs coupling to quarks after the quark mass 
has been factored out and QhiVnUf. is defined in Eq. 23 



Because of the dependence on the Higgs mass and the fact that the Higgs couplings 
to quarks goes only through the doublet component, the lightest Higgs doublet generally 
gives the dominant contribution. Note that the Higgs contribution is inversely propor- 
tional to the square of the sneutrino mass and is therefore expected to dominate at low 
masses since the Z contribution depends only weakly on the sneutrino mass through the 
effective mass, ^un- Furthermore the Higgs contribution is roughly the same for neutrons 
and protons. The quark coefficients of the nucleon were taken to be the default values 
in micrOMEGAS2.4 (with /^^ ^ = 0.033, 0.023, 0.26 and /^^ , = 0.042, 0018, 0.26) l69 . 
There can be large uncertainties in these coefficients, in particular recent lattice results 



point towards a smaller s-quark coefficient 70 , however the quark coefficients will have 



a significant impact only when the Higgs contribution is dominant, that is for sneutrinos 
below ^ 100 GeV. 

The total spin independent cross section on point-like nucleus is obtained after averag- 
ing over the sneutrino and anti-sneutrino. Note that the interference between the Z and H 
exchange diagrams have opposite signs for vN and i'*N . Here we assume equal numbers 
of sneutrino and anti-sneutrinos so that the total cross section is the sum of the Z and 
H contributions. To compute the direct detection cross section we use micrOMEGAs2.4 
and thus include all the contributions. To take into account the fact that the proton and 
neutron contribution are not necessarily equal, we compute the normalized cross-section 
on a point-like nucleus. 



SI _ y^JP ' V' ^)Jri) . „s 

where the average over and v\ is assumed implicitly. This cross-section can be directly 

XP 



compared with the limits on cr|p that is extracted from each experiment. 



5 Results 

The free parameters of the model are m^^, /x. Mi, M2, M(, Mk-, ^a, ^Zi-, Oe(,,oiz, as well as 
all masses and trilinear couplings of sleptons and squarks. To reduce the number of free 
parameters we fix those that do not belong to the sneutrino or neutralino/Higgs sectors. 
We thus fix the soft masses of sleptons and squarks to 2 TeV and we take Af = Q and 
At = ITeV. We furthermore assume Mi = M2/2 = M3/6 as dictated by universality at the 
GUT scale and Mk = 1 eV since this parameter mainly affects the B — B' mixing and is 
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not directly relevant for our study. The remaining free parameters are m^^, fi, Mi, M[, Ax 
and Mz2, oiz- We first consider specific choices of 6'Eg before letting it be a free parameter. 

5.1 The case of the UiX)-^ model 

As we have discussed above, we expect the relic abundance of the sneutrino dark matter 
to be generally too high. The only ways to bring the abundance within the range preferred 
by WMAP will be to enhance the annihilation through a resonance effect or make use of 
coannihilation. We will first describe the typical behaviour of Q/i^ as a function of the 
vr mass for fixed values of the free parameters. Once Qe^ is fixed the parameters that 
have a strong influence on dark matter arc Mz.^^, the mass of the new gauge boson, A\ 
that influence the Higgs spectrum, as well as \i, Mi that determine the region where the 
sneutrino is LSP through their influence on the neutralino mass. The parameter /i also 
influences the Higgs mixing matrix and therefore the coupling of the Higgses to sneutrinos. 
To limit the number of free parameters we will first fix Mi = M[ — 1 TeV and rather 
modify the property of the neutralino NLSP by varying /x. Furthermore allowing to have 
fi < Ml will cover the case of the higgsino NLSP which has much higher annihilation 
cross section than the bino and is therefore more likely to play an important role in 
coannihilation channels. 

5.1.1 A case study with Mz^ = 1.6 TeV 

We first consider the case where the new gauge boson is just above the LHC exclusion 
limit, that is Mz2 = 1-6 TeV and we fix = 1 TeV, az = —0.001 and Ax = 1.5 TeV. 
For this choice of parameters, the mass spectrum is such that /i2 is dominantly singlet 
and has roughly the same mass as the Z2 while the other heavy Higgses (/is, ^4°, /i"*") are 
nearly degenerate with a mass around 1.97TeV. The light Higgs has a mass of 136.6 GeV 
and is safely above the LEP hmit despite the small value of tan^ = 2.05. The lightest 
neutralino has a mass of 957 GeV and is a mixture of bino and higgsino. The sneutrino is 
therefore the DM candidate when its mass ranges from a few GeV's up to the mass of the 
lightest neutralino. The value of Qh'^ lies below the WMAP upper bound in two regions, 
the first when 60GeV < m£>^ < 69GeV the second when 734GeV < m^^ < 805GeV. 
In both regions the annihilation cross section is enhanced significantly by a resonance 
effect. In the first region the hi exchange is enhanced and the preferred annihilation 
channel is into bb pairs while in the second region the h2 and Z2 exchange are enhanced 
with a dominant contribution from h2- The preferred annihilation channels are into 
W'^W~ , ti, hihi, ZiZi pairs, the dominant decay modes of /i2- Note that the decays into 
gauge bosons proceed through the small doublet component of /i2- For larger masses of 
the sneutrino, coannihilation with the lightest neutralino can take place. Coannihilation 
processes involving pairs of neutralinos annihilating through the heavy pseudoscalar Higgs 
can decrease flh^ such that the WMAP upper bound is satisfied. This occurs when the 
NLSP-LSP mass difference drops below 60 GeV. 

This general behavior is somewhat influenced by our choice oi Ax, /J, and az since these 
parameters influence the masses and mixings of the Higgs doublets and the couplings of 
sneutrinos to Higgses. For example for Ax = 100 GeV, the mass of the heavy doublet 
(which is now /12) decreases to 592 GeV. The h2/A^ exchange contributes to sneutrino 
annihilation into top pairs thus leading to a decrease of the value of Vlh? as one approaches 
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Figure 3: Qh"^ as a function of the LSP mass for Mz^ = 1-6 TeV and (a) 
/i = 1 TeV, {AxiTeV),az) = (1.5, -IQ-^), (0.1, -10-=^), (2., -10"^) (b) = l.STeV, 
{AxiTeV),az) = (1.4, -10-=^), (1.4, -10"^) and /i = 0.7TeV, {Ax{TeV),az) = 
(0.6,-10^^). Full (dash) lines correspond to the region where the LSP is the RHSN 
(neutralino) . 



the Higgs resonance, the drop is not significant enough to bring the relic density below the 
WMAP upper bound, see Fig. [3^. As another example consider the case A\ = 2 TeV and 
az = 10""^. Here the hiVRVR coupling is suppressed, fi/i^ becomes very large and despite 
a resonance effect when m^^ ~ m/jj / 2 the WMAP upper bound can never be satisfied for 
a sneutrino hghter than lOOGeV, see Fig. |3^. 

The parameter induces corrections to the light Higgs mass as well as shifts in the 
Higgs mixing matrix. In particular increasing /i (and therefore A) to 1.5 TeV increases the 
singlet mixing in the light Higgs and thus the vrVrHi coupling. This makes annihilation 
processes through Higgs exchange more efficient, and gives rise to a new region where 
VLh'^ is below the WMAP upper bound when the sneutrino mass is just above the W pair 
threshold, see Fig. |3)d with Ax = 1.4 TeV. Note that in this case the lightest neutralino 
has a large bino component and coannihilation is not very efficient. Lowering fi to 700 
GeV, changes the nature of the lightest neutralino which becomes dominantly higgsino 
with its mass determined by /z. Thus the range of masses where the sneutrino is the 
LSP becomes narrower. In fact the singlet Higgs/Z2 pole annihilation region cannot be 
reached when the sneutrino is the LSP, the region compatible with WMAP is rather one 
where higgsino coannihilation dominates. In Fig. |3|d, one can see a significant drop in 
Qh"^ near the /12 resonance, this is however not sufficient to have efficient annihilation of 
the RHSN. 

The direct detection spin independent (SI) cross section receives contributions from 
both the light Higgs and Zi exchange. In general for the U{1)^ scenario, the cross- 
section is below the limits from XenonlOO ff]. As we have argued above, this is because 
the Zi contribution is directly proportional to sin^a^. This contribution nevertheless 
dominates for a mixing angle \az\ = 10~^ except for small masses. It gives a = 1.4 x 
10~^° pb for m^;^ > 200 GeV and is mostly independent of other input parameters, 
see Fig. |4^,b. Furthermore the Zi contribution is much larger for neutrons than for 
protons. We take that into account by computing the cross section for scattering on point- 
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Figure 4: crf^^ as a function of the LSP mass for the same parameters as Fig. 3 The 



XenonlOO exclusion curve is also displayed. 



like nucleons (here we use Xenon) normalized to one nucleon. For masses in the range 
50—100 GeV, the SI cross section is enhanced due to the light Higgs exchange contribution 
which increases at low sneutrino masses, although for this scenario the values are always 
below the experimental limit. The predictions are in the range a = 2 x 10^^*^ — 3 x 10^^° pb 
for a = — 10~^. The SI cross section drops rapidly with the mixing angle, for example for 
az = |10~^|, o'fj^AT — 10"^^ pb for m^^^ > 200 GeV, see Fig. For small values of the 
mixing angle, the Zi contribution is suppressed and the Higgs contribution can become 
dominant even for masses of a few hundred GeV's, this however corresponds to a low 
overall cross section, see for example Fig. |4)d. The spin-independent cross section can be 
further suppressed in the limit az = and for parameters for which the hi couplings to 
sneutrino is suppressed. 



5.1.2 Exploration of U{1)^ parameter space 

After having describe the general behavior of Qh"^ and of the direct detection rates for 
some choice of parameters, we now search for the region in parameter space that are 
compatible with 0.1018 < ilh'^ < 0.1228 corresponding to the WMAP 3cr range [tT 



as 



well as with the direct detection limit of XenonlOO [T] . We have further imposed the limits 
on the Higgs sector from LEP and the LHC and have taken into account the impact of 
the invisible decay mode of the Higgs. We have also used limits on sparticles from LEP as 
implemented in micrOMEGAs . Note that since we are only considering the case of heavy 
sfermions, these limits as well as the new LHC limits on sparticle masses do not come 



into play [59 ,60 . Finally we have also imposed the constraints from AM^^s- We have 
performed random scans with 5 x 10^ points. The range for the parameters are listed 
in Table |2} As before we here fix all sfermion masses to 2TeV and neglect all trilinear 
couplings with the exception of At = I TeV. Note that we impose > and therefore 
consider only positive values for Ax, see Eq. [15] 

The results of the scan as displayed in Fig. [5] show that the allowed values for rrij^^ 
cover a wide range from 55 GeV to the largest value probed. The allowed points in 
the plane — Mz^ are clustered in three regions around ^ 60 GeV, around 
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Parameter 


Range 




[0, 1500] GeV 


Mz, 


[1300, 3000] GeV 


fi 


[100, 2000] GeV 


Ax 


[0, 2000] GeV 


az 


[-0.003, 0.003] rad 


Ml, M[ 


[100, 2000] GeV 



Table 2: Range of the free parameters in the U{1)' models 



mp^ ^ Mz^ - S{+S') with 5 ^ 70 GeV{S' ^ 10 GeV), and m^^ ^ 90 - 100 GeV. The first 
two are characterized by the main annihilation mechanism near a resonance. The latter 
corresponds to annihilation through a Higgs exchange. As discussed above, this requires a 
large value for /x, see Fig. [SjD. The remaining allowed scenarios correspond to mi>^ ~ ^xi- 
In most cases the NLSP is a higgsino and m^^ ~ fi, see Fig. |5jD, then important contribu- 
tions from coannihilation processes involving neutral and charged higgsinos annihilating 
into fermion pairs are to be added to the sneutrino annihilation processes dominated by 
the channels into WW, ZZ through Higgs exchange. Because /i is constrained by the 
LEP limit on charginos, the RHSN mass is in this case heavier than 90 GeV. A few cases 
where the NLSP has an important bino component are also found, those correspond to 
the points above the line /x ~ m^;^ in Fig. |5}d. 




M, (TeV) (TeV) 



Figure 5: The allowed scenarios in the a) Mz2 vs m£>^ plane and b) fi vs m^^ plane. Points 
excluded by the AMd^s constraint are displayed in red. 

In this model, the cross section for LSP scattering on nuclei is dominated by either 
Higgs or Zi exchange and is mostly independent of Mz2 ■ The direct detection cross section 
can reach at most cif^^ = 2 x 10~^ pb when m/>^ ~ 50 — 60 GeV and the light Higgs 
exchange dominates while the maximal value is only cr?\ ^ 5 x 10~^° pb for sneutrinos 
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above 200 GeV when Zx exchange is dominant. These predictions are much below the 
exclusion hmits of XENONIOO [7]. Values below af^^v = 10"^^ pb can also be obtained, 
these occur when the Z\ contribution is suppressed by the small mixing angle and the light 
Higgs coupling to the LSP are suppressed as well. Cases where coannihilation processes 
dominate can also lead to small cross sections. 

Constraints from AM^ rule out some of the parameter space, in particular the case 
where > since it leads to values of tan /3 < 1. The allowed points in the tan /3 — 



plane are displayed in Fig. mb. As mentioned in Section 2.5, when tan/? is small the 



charged Higgs box diagram adds to the SM contribution and leads to value for AM^ that 
is too large. This constrain does not influence directly the range of predictions for the 
direct detection cross section, nor does it affect the range of allowed masses for the RHSN. 




\m a^frad) 



Figure 6: (a) cr?l^xe ^ function of the LSP mass for the allowed scenarios in the U{1)^ 
model, (b) Allowed scenarios in the tan/3- m^^^ plane. Points excluded by Aikf^ are 
displayed in red. 



5.2 The case of the U{l)r] model 

The properties of the RHSN dark matter are dependent on the choice of the U(l)' charges. 
To illustrate some of the differences we present results for = — .297r (f/(l)^) before 
discussing arbitrary values in the next section. We chose this model for its phenomeno- 
logical properties. For this value of 6'Eg, the coupling of the sneutrino to Z2 is enhanced 
as compared to the previous example with Q'^, = — -\/5/2a/3. Furthermore the coupling 
of the RHSN to hi is typically enhanced as compared to the model f/(l)^. Finally the 
vectorial couplings of Z' to neutrinos and d-quarks are non-zero. Therefore both terms in 



Eq. [24] can contribute to the SI cross section, this implies a direct detection cross section 
that is much larger than found previously. Nevertheless the vectorial couplings of Zi are 
far enough from their maximal value that it is possible to find scenarios that satisfy direct 
detection limits as will be demonstrated below. 
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5.2.1 A case study with Mz^ = 1.6 TeV 

To illustrate the behaviour of the relic density and the direct detection rate we choose 
= 1-6 TeV, az = 0.001, /i = 1.5 TeV and A\ = 0.5 TeV. For this choice of parameters 
rrihi = 116.9 GeV despite a small value of tan/3 = 1.2. The doublet Higgses are around 
1.3 TeV while is nearly degenerate with Z2. As for the the model U{1)^, the relic density 
satisfies the WMAP upper bound in the regions where annihilation into light Higgs or 
singlet Higgs/Z2 is enhanced by a resonance effect as well as in a region where annihilation 
via hi exchange is efficient without the benefit of a resonance enhancement. The latter 
region extends over a wide range of values for the LSP mass because of the large couplings 
of the sneutrinos to hi. The preferred annihilation channels are typically into WW, ZZ 
tt or hihi pairs as well as into bb for light sneutrinos. Note also that exchange of /12 can 
contribute significantly to sneutrino annihilation, see the small dip at m,>^ = 650 GeV in 
Fig. [7^. In this scenario the direct detection cross section is large (cf^^ ^ 4.5 x 10~^ pb 
for m,>^ > lOOGeV) and exceeds the XenonlOO bound for a sneutrino LSP below 550 
GeV, see Fig. [7]d. This is because both the Zi and Z2 contribute significantly to the cross 
section. Extending the region of validity of sneutrino DM over a larger mass range after 
considering direct detection limits thus requires decreasing the mixing angle az and/or 
increasing the mass of the Z2. For example taking az = 10~^, fi = Ax = l.STeV has 
the effect of decreasing the relic density - so that it is below the WMAP upper bound 
for sneutrino masses in the range 90-900 GeV - while also decreasing the SI cross section. 
Yet the LSP is still constrained to be m,>^ > 370 GeV. In fact the direct detection rate 
decreases by less than a factor of two, this is because the contribution of the Z2 exchange, 
the term proportional to Qy in Eq. 24, is not suppressed by sina^. This means that to 



have sneutrino DM with a mass of a few hundred GeV's also requires increasing the mass 
of the Z2. For example taking Mz2 = 2.5 TeV decreases the direct detection rate below 
the XenonlOO exclusion for all masses of the LSP. However for this choice of parameters, 
only light sneutrinos also satisfy the relic density constraint since the heavy sneutrino 
LSP cannot annihilate through the singlet/Z2 resonance. 

Finally we comment on the coannihilation region. We have fixed Mi = 1 TeV, so 
that the lightest neutralino has a mass ~ 1 TeV and is dominantly bino. For values of 
/i < 1 TeV, the lightest neutralino has a large higgsino component and coannihilation is 
efficient. For the benchmarks with fi = 1.5 TeV the coannihilation processes involving 
binos are not efficient enough to reduce Qh"^ to a value compatible with WMAP. 



5.2.2 Exploration of t/(l),, parameter space 

We also explore the parameter space of the model, varying the parameters in the range 
shown in Table [2j As before, we impose in addition to limits on fi/i^ and on the SI cross 
section, the lower limit on the Higgs and Z2 mass, and the limit on AM^^s- The results 
are displayed in Fig. [8] in the Mz^ — rrij^^ plane, as well as in the fi — rrij^^ plane. As we 
have discussed above, because Higgs annihilation is efficient enough, the sneutrinos are not 
confined to the hi, Z2 resonance region. Sneutrinos from lOOGeV to Mz^l'^ can satisfy the 
relic density constraint, either through pair annihilation or through coannihilation (the 
region where m^^^ ~ /^). In the former case large values of are preferred to have large 
enough couplings of the sneutrino LSP to hi. The direct detection cross section is large 
with a maximum value near cr|^^ 10~^ pb for the whole range of LSP masses. Because 
the experimental limits on the direct detection rate are more stringent for masses around 
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Figure 7: flh"^ as a function of the LSP mass for a) = 1-6 TeV, /i = 1.5 TeV, 

Ax = 0.5 TeV, az = 10"=^ and Mz^ = 1.6,2.5 TeV, fx = 1.8 TeV, Ax = 1.8 TeV, 
az = 10~^ . The dotted hues correspond to the neutrahno LSP. b) crf^^ as a function of 
the LSP mass for the same choice of parameters. 



80 GeV, hght sneutrinos are severely constrained especially when the Z2 mass is lighter 
than 2TeV. The lower bound on the cross section for < 3 TeV is = 2. x 10~^ pb 
. The typically large direct detection rate is the main signature of this scenario. To 
completely probe this scenario up to DM masses of 0.3(1.5) TeV requires one (two) orders 
of magnitude improvement over the current limits. Several of the scenarios are also 
constrained by AM^, in particular when az < which implies tan/3 < 1. This does not 
affect the range of values of sneutrino masses compatible with all constraints. 




(TeV) M, (TeV) 



Figure 8: The allowed ?7(1),, scenarios in the a) vs m,>^ plane and b) /i vs rrij^^ plane. 
In red, points excluded by the AM^^g constraints, in green those excluded by XENONIOO 
and in yellow those excluded by both constraints. 
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Figure 9: a) o"f^xe a function of the LSP mass for the allowed scenarios in the f/(l)r) 
model, b) Allowed scenarios in the tan (3- az plane, same color code as Fig. [sj 

5.3 A global scan of the parameter space 

Having illustrated the properties of the sneutrino dark matter in specific models we will 
next explore the complete parameter space of the model by keeping all parameters for 
the neutralino, sfermion and gauge boson sector free while imposing the constraints from 
WMAP, direct detection as well as on Higgs masses. 

The values of where the sneutrino is a good dark matter candidate are restricted. 
First for 6'£;g ~ the value of Vs becomes very large especially when is large. This 
induces large negative corrections to sfermion masses and lead in particular to a charged 
LSP. For example, for soft terms at 2TeV the values —0.2 < 9e^ < 0.05 are excluded when 
= l.TeV. Second the direct detection cross section is often too large when \OEfi\ < 0.5. 
This is due mainly to the contribution of the Z2 exchange to the direct detection cross 
section that is proportional to cosOeq- To illustrate this we display the variation of the 



direct detection cross section as a function of OEa for different values of M^^ in Fig. 10 
We choose m,>^ = Mz^fi so as to guarantee that the WMAP upper bound is satisfied and 
fix /i = M^2/2 + 0.5 TeV in each case to ensure that the sneutrino is the LSP. We have 
also fixed Mj = 2 TeV. The direct detection bound are easily satisfied near Oe^ = ±7r/2 



because of the suppressed contribution of the Z vectorial coupling, see Fig. [TOj The mixing 
angle az has only a moderate impact on the direct detection rate, while increasing the 
mass of the Mz^ reduces the SI cross section, except when OEa = 0^ as we have seen in 



Section 5.1 Note that for < there is a dip in the cross section near Oeh = 7r/2, this is 



because of a cancellation between the y and y' contributions in Eq. 24 in the cross section 
on neutrons. 

For the random scans we choose the same range as in Table [2] with in addition — 7r/2 < 
(^Ee < 7r/2. As before, we have applied the constraints on Mz2 from the LHC and on az- 
We have imposed the latest constraint from XENONIOO and from AMg a posteriori 
to better illustrate the impact of these constraints. The successful scenarios are displayed 



in Fig. 11 in the Mz2 — ^e^ plane. The white region at small values of Oeq have a charged 
fermion LSP. Many scenarios with —ir/A < 9eq < have too large a direct detection 
rate as was the case for the model f/(l),,. The value of the direct detection cross section 
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span several orders of magnitude from less than 10^ pb to 2 x 10^ pb. In particular 
for sneutrino masses around 100 GeV there are many models which exceed the direct 
detection limit. These are scenarios with Oeq < ^ where the sneutrino coupling to light 
Higgs allows efficient annihilation even away from resonance as discussed in the case of 
U{l)fj. The enhanced couplings to the Higgs as well as the coupling to the lead to 
a large direct detection rate. Models with 6'^;^ > have only light sneutrinos in that 
mass range when coannihilation plays a role and therefore tend to have much lower direct 
detection rate. In fact even for heavier masses models with Oe^ > predict smaller direct 
detection rates as discussed above. 




Figure 11: a) Allowed scenarios in the plane, the models above the XENONIOO 

bound are in green (light grey), b) crf^jy as a function of 

The processes that can contribute to flh"^ were discussed both in the context of U (1)^ 
and U{l)rj models. In the general case we find similar results. We find a predominance 
of annihilation near a hi or singlet Higgs/Z2 resonance as well as annihilation into gauge 
boson pairs through hi exchange. The latter being confined to sneutrino masses just 
above the W pair threshold when Oe^ > 0. These regions have a high density of points 
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in the plane Mz^ — '^vr iii Fig. 12 For 6eq > ^ the only other allowed scenarios have 
displayed in Fig. U2p. These are dominated by Higgsino coannihilation. For 
(^E(i < 0, sneutrinos of masses above 100 GeV can also annihilate efficiently through hi 
exchange provided their coupling to hi is large enough - this requires large values of fi. In 
both figures we have imposed the AMd^s constraint although we do not display explicitly 
its impact. As discussed in the previous section this constraint plays a role for small 
values of tan P for any values of the LSP mass. 




6 Conclusion 

The right-handed sneutrino is a viable thermal dark matter candidates in U(l) extensions 
of the MSSM. The allowed parameter space depends strongly on the value of the new 
Z2 vector boson mass. Sneutrino annihilation is typically dominated by resonance anni- 
hilation, with in particular the dominant contribution from the Higgs sector rather than 
from the new gauge boson Z2. For the light Higgs this requires fine tuning of the masses 
while for the heavy singlet Higgs the mass difference has to be within 15% (of mhj2). 
As in other supersymmetric models, coannihilation processes can also be important. For 
simplicity we have only discussed the case of one sneutrino dark matter candidate, in the 
case of complete degeneracy of three sneutrinos the relic density increases by a factor 3. 
Indeed the annihilation cross section is the same for all sneutrino flavours since there is 
never a significant contribution from annihilation into neutrino pairs, thus the increase 
in the number of channels is compensated by the increase in the degrees of freedom lead- 
ing to a smaller effective annihilation cross section. This will imply a narrower range of 
masses for the LSP near Mh./2. 

The direct detection limit is very stringent for a whole class of models unless the mass 
of the Z2 is above 2TeV. These scenarios with ^_Bg < will be best tested with SI detectors 
with improved sensitivities. In particular with a factor of 2 better sensitivity the whole 
region with intermediate masses for the sneutrinos - when annihilation into W pairs is 
dominant - will be probed. The scenarios with ^^jg > are more challenging to probe 
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via direct detection. Another signature of the RHSN model is that of a new gauge boson 
at the LHC In fact the recent results from searches for gauge bosons at the LHC 
have had a significant impact on the parameter space of the model and pushed one of the 
favored region for the sneutrino LSP to be near SOOGeV. A Z2 up to several TeV's will be 
probed when the energy of the LHC is increased to 14TeV. A negative result from such 
a search will imply a much reduced rate for direct detection searches as well. Indirect 
searches could also provide a good probe of the RHSN model, these will be investigated 
separately. 

In this analysis we have assumed that Dirac neutrinos were very light (sub-eV range). 
Such right-handed Dirac neutrinos can be produced by Z' interactions prior to big bang 
nucleosynthesis leading to a faster expansion rate of the universe and to too much He^. 
The resulting constraints on the mass of the new gauge bosons in the model were 
analyzed in 72 and compatibility with BBN resulted into lower limits on the Z' mass 



in the multi-TeV range assuming the effective number of neutrinos was increased by 0.3. 
However since then new data from WMAP have been released which implied a higher 
than expected value for the abundance of He^ and for the effective number of relativistic 
neutrino species, more precisely ^ejf 



73 , 74 . This new value would 



4.34 + 0.86 -0.88 

have the effect of relaxing the lower limit on the Z' mass to a value comparable to LHC 
bounds. 
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A The masses in the Higgs sector 



The tree-level CP-even Higgs mass matrix in the basis S") reads 53 
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The dominant radiative corrections due to top and stop are 



T s 




{MlU = k 



[m 



2\2 



im~ — m~ Y 



2V2 



where k = 3/(47r)^, ht is the top Yukawa couphng and 



g = 2 

J^ = log 
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2 2 

m~ — rrix 

tl t2 

2 2 
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log 



and 



ml = h'^ji {ji - Atti3) 



ml = h^tAt At - 



(29) 



(30) 



(31) 



The DR renormahzation scale, A is fixed to 1 TeV. The CP-odd mass tree-level matrix 
reads 



(32) 
(33) 
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^,JG {1,2,3} 



(34) 



where vi = Vd, V2 = f«, fs = f s- 

The charged Higgs mass including radiative corrections from the top and stop sector. 
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